
myps,'./emaf_paper/plots/hyperion_ffore_midplane.eps',/bw,xsize=3.1,ysize=2.5

;; Load in the ffore cube
cube = readfits('./local/MW_model_ffore.fits',hd)

plane = median(cube[*,19:20,*],DIM=2,/even)

sxaddpar,hd,'NAXIS2',sxpar(hd,'NAXIS3')
sxdelpar,hd,'NAXIS3'
sxaddpar,hd,'NAXIS',2

sxaddpar,hd,'CD2_2',sxpar(hd,'CDELT3')/1.d3
sxaddpar,hd,'CRPIX2',sxpar(hd,'CRPIX3')/1.d3
sxaddpar,hd,'CRVAL2',sxpar(hd,'CRVAL3')/1.d3

cs = 0.7

pr = [0,1]
plotmap,plane,hd,ct=0,range=pr,axsel=9,charsize=cs,$ ;xmargin=[10,10],$
        XC=xc,YC=yc,xmargin=[8,2],ymargin=[3,0]
cgcontour,plane,xc,yc,levels=(findgen(10))/10.,/over,c_thick=1.5,$
          c_colors=[replicate('white',7),replicate('opposite',3)],$
          c_charsize=cs*0.75

;; Edge of the survey area
vline,7.5,thick=3,linestyle=3
cgloadct,34,/silent

;; Tangent distance
ell = findgen(651)/10.
dtan = 8.5*cos(ell*!dtor)
cgOplot,ell,dtan,color='Opposite',thick=6,linestyle=5

boxedge = 15/cos( (90.-ell) * !dtor )
cgOplot,ell,boxedge,color='Black',thick=1.5

myps,/done



END
